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SUMMARY 



The author derives equat ions containing only directly 
measurable constants for the quantities involved in the 
hydrodynamic theory of detonation. The stable detonation 
speed, D, is revealed as having the lowest possible value 
in the case of positive material velocity, by finding the 
minimum of the Du curve (u denotes the speed of the gases 
of combustion). A study of the conditions of energy and 
impulse in freely suspended detonating systems leads to 
the disclosure of a rarefaction front traveling at a lower 
speed behind the detonation front; its velocity is com- 
puted. The latent energy of the explo si ve . pas se s into the 
steadily growing detonation zone - the region between the 
detonation front and the rarefaction front. The conclu- 
sions lead to a new definition of the concept of shatter- 
ing power. The calculations are based on the behavior of 
trinitrotoluene. 

INTRODUCTION 



The chemical change in an explosive substance, such 
as the propelling charge in a firearm, generally takes 
place at a linear combustion speed in fractions of meters 
per second, whereby this rate at which the conversion zone 
penetrates the explosive substance is a function of the 
gas pressure. 

Berthelot (reference l), in 1881, discovered the pe- 
culiar fact that under suitable conditions, substantially 
higher conversion speeds are feasible. In his experiments 
with explosible gas mixtures in long tubes, he ascertained 
propagation speeds in the conversion zone of several kilo- 
meters per second, whereby this speed remained constant 
over the entire conversion distance and even remained un- 
affected by the initial gas pressure. 

*"Beitrag zur l( hydrodynamischen Detonationstheorie . " Zeit- 
schrift fur technische Physik, vol. 19, no. 9, 19S8, 
pp. 271-283, 
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Berthelot and Dixon essayed a molecular kinetic ex- 
planation of the detonation processes and, in fact, achieved 
with many gases, practical values for the detonation speed. 

The purely thermodynamic and aerodynamic treatment of 
the detonation process which disregards special molecular 
concepts, is attributable to Schuster (reference 2), and 
Chapman (reference 3), who utilized the findings obtained 
by Riemann (reference 4), and Hugoniot (reference 5) in 
the treatment of the linear compression shock in gases. 

The extension of these macroscopic theories to in- 
clude solid explosive substances, is due to Jouguet (refer- 
ence 6) and R. Becker (references 7 and 8), the latter 
dealing largely with the cause of detonation stability. 
He proves, with the help of the thermodynamic entropy laws, 
the validity of Chapman's term for detonation speed. 
Schmidt's more recent studies of the detonation process of 
solid explosive substances (reference 9), containing a 
wealth of experimental material, are particularly well- 
suited to prove the validity of the hydrodynamic theory of 
detonation. 

However, there is a gap in the orthodox theory. One 
obtains for the combustion gases behind the detonation 
front, a speed in the direction of the advancing front. 
The consequences of the existence of such a speed on the 
total system, are not considered for the reason that only 
the measurable quantity of this speed of the combustion 
gases is of influence for the detonation opoed (reference 
10 ). (Schweikert refers to this difficulty, but is of the 
opinion that it would void the whole hydrodynamic theory 
of detonation. ) 

To illustrate the defect of the present theory, the 
derivation of Chapman and Becker is briefly reviewed. 



CHAPMAN AND BECKER EQUATIONS 

A tube of 1 square centimeter section contains an ex- 
plosible mixture under pressure p 0 , density 'p , specif- 
ic volume v , and temperature T 0 (fig. l). 

The front of the conversion moves toward the right in- 
to the explosive at the empirically constant detonation 
speed D. In zone JP the chemical change takes place, 
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the product of which, the gases of .combustion, also moves 
toward the right at speed u. 

For an observer, at rest relative to the combustion 
gases, hence also moving rightward at speed u, the pres- 
sure Pj , density p L , specific volume v , and temper- 
ature T x , prevail- Certain relations must exist "between 

these macroscopic quantities which result from the general 
physical laws governing the preservation of energy, im- 
pulse, and mass. 

Visualize an observer in the conversion zone S 1 moving 
at detonation speed toward the right. To him the process 
of conversion must appear stationary if D is constant. 
The explosive moves at speed D into zone P from the 
ri^ht, while the coin "bust ion gases leave the zone at speed 
D - u toward the left. 

The preservation of the mass gives: 

D p 0 = (!) - u) p x (1) 

or 

L v x = ( D - u ) v 0 

The momentum of the explosive substance moving into 
the zone per second from the right is 

B Q = p Q D D 

The momentum of the combustion gases leaving the 
zone per second is 

3 X = p x (D - u) (D - u) 

The difference B - B -of both momentums must be 

o 1 

equal to the effective pressure difference. Hence, 

B o - ?i = P x - Po 



and, in connection with equation (l), the momentum equa- 

} (2) 



tion 

p - t> = p D u 
*i -o r o 



(P, - Pn> V n = D U 



s 
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From the right, four different, forms of energy enter 
the zone per second: 

1. The latent chemical energy. If this amounts to 

U per unit mass, the leftward moving energy per second is 

Po D U 

2. The kinetic energy 

3. The heat contents of the undecomposed explosive 

gas 

P 0 D c v T o 

if c y denotes the specific heat at constant volume per 
unit mass. 

4. The external energy 

Po D 

Leaving the conversion zone per second, toward the 
lef t , are : 



1. The kinetic energy of the combustion gases: 
| Pl (D - u) ( E - u) 2 



2. The heat content of the combustion gases 

p 1 (D - u) c v T x 

3. The external energy 

Pj (D - u) 

The law of the preservation of energy therefore gives 
the next equation 

P 0 D U + \ p Q D D 2 + p Q D c T T Q + p Q D 

= | Pl (D - u) (D - u) 2 \ (3) 

+ p i (D - u) c v T x + p x (D - u) 



i 
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or, in conjunction with equations (l) and (2) ^ 

| (P, + P o ) (v Q - v x ) = c v (T i - T q ) - U "~ (4) 
whereby 

and 

D = v 0 /HZH I 
0 J v o - v i 

as is directly obtainable from equations (l) and (?). 

Now the substitution of p, and v, or p 0 and v 0 
for T 1 and T 0 in equation (4) through the phase equa- 
tions, gives the relation between p x and Vj, with the 

constants of the original gas as parameters. 

y 

This relation between r> 1 and v 1 is the so-called 
"dynamic adiabatic" or Eugoniot curve, named after its in- 
ventor (fig. 2). 

Accordingly, if in a gas of original state, p 0 and 
v 0 , a conversion zone travels at constant speed - saying 
nothing at present about the magnitude of this speed - the 
state of the combustion gases is represented by some point 
on the Hugoniot curve. According to equation (6) every 
point has a certain speed of propagation D e.nd., accord- 
ing to equation (5), a certain speed of combustion gases. 

For all points above G-, there is a speed of combus- 
tion gases in the same direction as that of the propaga- 
tion of the conversion zone. This conversion is associated 
with a marked pressure rise and passes very quickly. If 
the conversion of the explosive has the characteristics of 
this zone, it denotes detonation. 



(5) 



(6) 



If the explosive mixture burns - i.e., if the speed 
of transformation is low, whereby the combustion gas pres- 
sure, p decreases relative to p 0 and the speed of the 
combustion gases is contrary to the direction of combus- 
tion - then the stage of the combustion gases is represent- 
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ed "by the curve "below E. The zone G- - E cannot he prac- 
tically realized. 

So, even though it is true that at detonation the stage 
of the combustion gases, and hence the magnitude of detona- 
tion speed, is characterized hy a point on the Hugoniot 
curve above G-, this still does not determine which point 
defines the actually occurring detonation speed, 

The only stable point on the curve, according to 
Chapman and Becker, is the contact point J of the tangent 
from A. 

For angle cp of the tangent , it is 



tan cp =, 




or, according to equation (6), 



33 = v 0 J tan cp (?) 



Since tan cp is computable from the dynamic adiabatic, 
expressions for D and u containing only macroscopic, 
thermodynamic equatioms can he adduced. 

The se are : 



■j k 

The constants 

k, ratio of specific heat 
E, gas constant 
M, molecular weight 
apply to the temperature T x of the combustion gases. 
Forming the difference a = 33 - u gives 
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. ..... a = y k — m~ 

But this expression presents, the magnitude of the velocity 
of sound in the combustion gases. 

Accordingly, it may "be stated: The detonation speed 
is given through the sum of the speed of combustion gases 
and sonic velocity in the combustion gases. 

The detonation temperature necessary for computing D 
is obtained from the specific heat and the heat content of 
the explosive whereby the temperature rise due to . the den- 
sity increase on passage through the wave front is still 
accounted for through the dynamic adiabatic curve. The 
detonation speeds computed this way are in good agreement 
with practical experience, hence affording conclusions re- 
garding the fundamental proof of the concept. 

^ The hydrodynamic derivation of formulas for solid and 
liquid explosives follows the same method, save that the 
specific volume of the gas molecules must be allowed for, 
as is customary in internal ballistics, by the introduc- 
tion of Abel's equation, a modification of van dor Waal's 
form, 

Quantitative calculation of. the detonation speed of 
solid explosives being out of the question, on account 
o"f ignorance of the covolumo of gases at high pressures 
and temperatures, the measured detonation speed is used 
instead for calculating the covolumes. 

In any case, the derived expressions indicate an ex- 
perimentally confirmed relationship with density p Q , 
which does not obtain for the gas formulas (3 ) and (9). 

DERIVATION OP HEW EXPRESSIONS 



7 

(10) 



The first task is the derivation of expressions for 
detonation speed D, and material speed u, which are 
■free from the concept of detonation, temperature and con- 
tain solely constants directly amenable to measurement. 
The derivation is to be carried out direct for solid and 
liquid explosives. By disregarding the covolume, the equa- 
tions become valid for gaseous explosives. 
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In the following, various concepts of internal "ballis- 
tics are employed.* 

The combustion of a quantity of explosive of weight 
L in a closed "bomb of constant volume V creates a ter- 
minal pressure p-g, called explosion pressure. 

For the covolume of gases per unit weight a, Abel's 
equation gives 

or, for the specific weight of the gases — = p 0 , 



„ _ * Po 

P B - f 1 - cx P( 



(11) 



f and a can be determined from pressure measurements in 
the bomb by means of equation (ll). 

The phase equation 

' N R T E rt 
PE (1 - a p Q ) = -g — P 0 

discloses that the experimentally defined constant f has 
the value 

* . 5^ («) 

with Tj;, explosion temperature. Tjj is the temperature 
of the. explosion gases existing when combustion takes place 
without external energy input; i.e., for instance, in the 
closed bomb. 

The gas under explosion pressure, is to be adiabat ical- 
ly released up to zenro pressure by means of a plunger. With 
k as exponent of the adiabatic release, the energy per- 
formed on the piston is 

*H. Langweiler, Z.f.d. ges. Schiess- u. Sprengstof fwesen : 
Development and experimental check of closed mathematical 
expressions for internal ballistic constants, particularly 
for the maximum gas pressure and the muzzle velocity. (To 
be published. )• 
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p E (V - a L) 

A = ~ - (13) 

k - 1 

which follows directly by integration of the adiabatic 
curve, according to Poisson. The energy capacity A of 
the gas in the form of equatfon (13) is called the poten- 
tial gas energy. But this energy A must be equal to 
the latent energy contained in the quantity of explosive 
L before combustion. If this amount per unit weight is 
U, we have 

p™ (V - a L) 

A = * I U 

k - 1 

or, with observance of equation (ll): 

f = (k - 1) U 

k = 1 + ^ (14) 

Constants f and U each have the dimension of an 
energy (m kg) per unit weight (kg), hence the dimension 
of a length (m). The exponent k i s an average value; 
but differs very little from that ( ~ 1 percent) computed 
from the specific heat, a,s. will be seen from the example 
later on. 

Other than density p Q of the explosive and co volume 

a x of the detonation gases, the expression for the deto- 
nation speed is to contain only the measurable constants 
U an d f . 

Figure 3 is another sketch of the detonation process. 
For the sake of clearness the two sides of the conversion 
zone carry pistons which at pressure p Q and p x act on 
the static explosive or the moving combustion gases. The 
conception of such pressure forces is necessary as other- 
wise no stationary phase would be conceivable. 

As the conversion zone moves toward the right for a 
distance D, the combvistion gases move for a distance 
u, in the direction of the advancing zone. Since the 
pressure p is to be constant within the combustion gases, 
the piston A likewise must travel a dists.nce u, toward 
the right. 
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The mass, impulse, and energy -equat ions for an ob- 
server moving with the detonation front, can he written 
in known manner. The pressure p 0 "being very small rela- 
tive to the detonation pressure p x (p G = 0 for solid ex- 
plosives), its effect can "be disregarded. 

Mass: 

D p Q = (D - u) p i (15) 

Impulse : 

p = ~ D u (16) 
l S 

(g is acceleration of gravity) 



Energy: 



D p ft U + - 1 - p n D ID 2 
0 2g 0 



p x (D-u) (l-a 1 p x ) 



k - 1 



(17) 



+ i p (D-u) (D-u) 2 + p (D-u) 
tig i l 



The left-hand side of the last equation gives the la- 
tent and kinetic energy of the explosive movement per sec- 
ond into the zone from the right. 

Behind zone F , leftward, a gas space of magnitude 
D-u per second is created, wherein the pressure p x 
prevails "because the piston A, as seen from the detona- 
tion front, moves a distance D-u to the left per sec- 
ond. How the potential gas energy of this newly created 
space is (similar to equation (13)): 

rs (D - u) (1 - a p ) 
- i l l 

k - 1 



The factor (1 - a x p x ) accounts for the reduction in gas 
space due to the co volume of the gases. 

The total gas space moves at a speed (D - u) toward 
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•the left. Its additional kinetic energy, accordingly, is: 

^ P i (D - u) (D - u) 8 

As seen "by the observer, the work p x (D - u) is per- 
formed on the plunger. The three forms of energy appearing 
per second on the left side of the detonation front must, 
together, equal the energy per second entering from the 
right. The result is equation (17). 

Quantities p x and P x may "be substituted herein 

through (15) and (16). With consideration of U (k - l) = 
f , we have : 

(k - l) D 2 + f g = fc-g-l (D - u) S 

+ (D - u) u (l - p Q + (k-l)u(D-u) 

and through direct calculation 



D = 



1 (l-fL + k_±_l) 

1 - a x p 0 \ u 2 J 



or abbreviated: 

Tl = 1 - a x p Q 
For negligible covolume T\ = 1. In 



n \ u 



K + k_±_l 
2 



(18) 



(19) 



(20) 



the detonation speed D is shown as function of the sole 
variable u, namely, the material speed. 

Since f has the dimension of length and g is the 
acceleration, the validity of the dimension. is proved. 

Figure 4 illustrates the function (20). The existence 
of a minimum is immediately apparent. The correlated pair 
of values is obtained by .differentiation of (20): 
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12 
du 



- 0 - i f- ^ + k_+_l N ) 



Hence 



u 



and, with (20) 



D = 



or 



k + 1 



f g 



1 



(21) 



</2(k + 1) f 



(22) 



D = £ (k + 1) u 



(23) 



Equations (21 ) and (22) are equivalent to Chapman's 
equations (8) and (9) as closer examination shows. Aside 
from the effect T| of the covolume, they merely contain 
the constant f known from bomb tests and the latent en- 
ergy U, "because k itself is = 1 + f/U. 



From the agreement of the 
detonation speed, according to 
stable from the Hugoniot curve, 
can be drawn: 



value of the lowest possible 
(20), with that derived as 
the following conclusion 



The constant detonation speed ultimately attained in 
a detonation process, is the lowest conversion speed pos- 
sible, according to the laws of physics, provided that a 
sufficiently strong initial ignition creates a combustion- 
gas speed in the direction of the advancing zone of trans- 
formation. 



The magnitude of the sonic velocity in the combustion 
gases is determined as follows: If E is the modulus of 
elasticity of the medium, the sonic velocity a in any 
medium follows the law: 



/" "E_g 
J P 



The modulus of elasticity for ideal gases is E = k p ob- 
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tained "fay. differentiation of the adiabatic curve 



.p 

But with 



p ^— J ts constant 



k 

p "* a i) = constant 



as. adiabatic curve of the combustion gases, its modulus of 
elasticity "becomes: 

p i k 
E = — 



1 " a i Pi 

Accordingly, the velocity of sound in these gases is 



a 



V P i (1 - a 1 p 1 ) 



and, after replacing p and p "by eauations (15) and 
(16) ~ -11 



/ k u (33 - u) 

J 1 - a i Po - 



D - u 

and, with the abbre viation T| = 1 - <x 1 p Q : 

a = i (k + 1) u - u (25) 

On comparison with equation (23), it is seen that the 
velocity of sound is equal to the difference between deto- 
nation speed 33 and combustion-gas speed u, that is: 

a = D - u! (26) 



The result is in accord with Chapman and Becker's findings. 
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At sonic velocity itself, the -equation is: 



or 



/k + 1 




k + 1 



f « 



(27) 



The detonation pressure can he computed according to 



equat ion (16 ) : 



v Po „ 

Pi = TT D u 
1 g 



The introduction of the values (21) and (22) gives: 



T) 2 f P o 



(28) 



Expressing the gas pressure of the explosives result* 
ing in the original volume under the slow combustion of 
the explosive, again as explosion pressure p s , and the 

covolume for this pressure with a E , Abel's equation 

f 'pQ 

P E 1 - ajj Po 

is applicable. 

From equation (28) then follows: 



2 Pe T - 



cc P _ 
i r o 



(29) 



and with it, the important result that the detonation pres- 
sure is twice as great as the explosion pressure, disregard- 
ing for the time being the effect of the difference in co- 
volume. For gases, it is exactly: 



■d=2 
•* l 



(30) 
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For the combust ion- £&s density, it is, according to 
equation (15), after insertion of . the D and u value: 



P, = 



k + 1 



k + 1 - T| 



(31) 



The detonation temperature T x can he computed in 

the following manner. The combustion gases follow the 
phase equation: 



But, according to equation (12), it is: 



and consequently, 



Hence , 



T, = 



R _ _f_ 

M T E 



T 



E 



(32) 



P x and p may he expressed with D and u: 

*E 



p D u /D - u \ 
Ti a _0__ r ^ 



Substituting equations (21) and (22) for D ; and u, 
we have the simple relation: 



m 2_k__ 

1 " k + 1 



E 



(33) 



The detonation temperature T x is, independent of 
the co volume of the combustion gases, greater by the fac- 
tor 'T~™T than the explosion temperature T-g of the ex- 

It T 1 

plo sive. 
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The aim of the first problem, the derivation of the 
quantities incident to the detonation of any explosive, 
containing - besides covolume - only the constants f and 
U, is herewith reached. Obviously, equation (22) is ap- 
plicable for the determination of covolume at of the com- 
bustion gases, in case the detonation speed D is meas- 
ured. 

THEORETICAL PROOF OF TEE EXISTENCE OP 

■- 

AN INCREASING DETONATION HEAD 



The second problem is the study of the consequences 
of a stationary detonation process illustrated in figure 
3. The assumption of a plunger A acting externally at 
magnitude of detonation pressure P i and moving at com- 
bustion-gas speed u in direction of the moving detona- 
tion front is, as already stated, necessary for the expla- 
nation of the stationary processes in the transformation 
zone E. 

At the initiation of a detonation process, the plunger 
A is represented by the nascent gases of the initiator. 
The energy of these gases, however, would only suffice to 
drive the detonation front a very short way into the ex- 
plosive. Once the action of the initiator gases has stopped, 
the cause for the positive speed of the combustion gases 
must naturally cease and the detonation come to an end. 

Rienann stated that a shock wave in gases could be 
maintained permanently only through external forces. But 
since experiments indicate that a detonation, once initi- 
ated, travels with ^reat constancy over any long distance, 
the explanation of the detonation process is in need of a 
supplement. This necessity becomes especially evident in a 
simple impulse consideration comprising a whole detonating 
system. 

Visualize a long column of explosive freely suspended 
horizontally. Ignition is to take place at one end. Af- 
ter the detonation wave has reached the other end, the 
then-formed gas of pressure p x and density P i would 

have to move as a whole at material velocity u in the 
direction of the detonation front. Such a motion is, of 
course, impossible, because no external forces which could 
cause a total impulse, act on the freely floating explosive. 
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The effect of the initial charge can he disregarded, 
since the column may "be of any length. 

An energy "balance comprising the whole system, also 
shows the incompleteness^ of the initially employed presen- 
tation; for the total energy of the combustion gases (po- 
tential and kinetic) of a freely suspended detonating col- 
umn of explosive is, according to calculation, greater 
than the latent energy originally existing in the explosive 

In the calculation this difference is covered "by the 
external energy of the plunger, hut in the absence of ex- 
ternal energies, as is the case for a column of explosive 
left to itself, no constant detonation speed would proba- 
bly exist. 

The difficulties introduced with the impulse and en- 
ergy balance of explosive systems left to themselves, must 
be voided by assuming a second wave front in the combus- 
tion gases. The speed w of this front relative to the 
quiescent explosive, must -be greater than the speed of the 
gases u, but smaller than the detonation speed. On pass- 
ing through this front the detonation pressure itself must 
drop to a value below that of the explosion pressure. 

Figure 5 illustrates the process of detonation of an . 
explosive substance on the assumption of a rarefaction 
front advancing at lower speed behind the detonation front. 
This rarefaction front moves at subsonic velocity relative 
to the combustion gases; that is, it is debatable, while a 
rarefaction shock of supersonic velocity would be impossi- 
ble. 

A - B is the explosive column made to detonate. Ig- 
nition takes place at A, where a solid abutment rests rel 
ative to the original explosive and which is therefore un- 
able, in contrast to a moving plunger, to perform work on 
the total system.. 

The significance of this "damming up" is referred to 
elsewhere. For the present, it may be stated that the 
omission of A modifies the speed of the rarefaction front 
but not the detonation speed. 

Fj_ is the detonation front wherein the pressure 
rises from p Q to the value p x . F 2 indicates the rare- 
faction front, within which p drops to p and the ma- 
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terial velocity changes from value . u to zero. Nothing 
is said about the form of the rarefaction front other than 
assuming that its extent is small compared to the distance 
P x F a . She" region "between the two fronts is called deto- 
nation head, and the distance F, F , the head length. 

Figure 5 indicates how, as a result of the difference 
in speed of the fronts, the length of the head continues 
to grow with time and, hence, with the detonation distance. 
Since the abutment A is fixed, the gas behind the rare- 
faction front F 3 must he at rest; i.e., its material ve- 
locity must become zero in case the process in the present- 
ed form is to be stationary. Suppose the values for this 
new stage of the combustion gases are: p 2 , p 2 , a s . 

At present the combustion-gas speed u is confined 
within the head F x F a . Before front Fj and behind 
front F 2 the material is quiescent. As the head becomes 
longer in unit time the momentum of the gases within the 
fronts increases with respect to time. The cause for this 
growth in. momentum lies in the pressure difference p - 

p Q existing at either side of the head. The pressure p g 

must be lower than the explosion pressure p-g , for the 

energy increment of the detonation head is necessarily at 
the expense of combustion-gas energy of stage 2. 

Assuming the nonexistence of the abutment at A, the 
inertia resistance of the gases in proximity of A is to 
be substituted. Then the gases behind the rarefaction 
front are subject to leftward acceleration - i.e., oppo- 
site to the gas speed within the head - so that the impulse 
.law is always complied with. Such leftward motion of the 
combustion gases in stage 2 is followed by a drop in pressure 
p relative to the pressure with fixed abutment. 

But a drop in p g calls for a decrease in temporal 
rise of gas momentum within the head; i.e., the head length 
would increase less rapidly in time or, ultimately, the 
speed of the rarefaction front would more closely approach 
the speed of the detonation front. This fact is of impor- 
tance for the effect of the explosive, as will be shown 
later on. 

In the following, the quantitative calculation of the 
speed of the rarefaction front with fixed abutment A is 
carried out. It could also be made for absent abutment, 



N.A.C.A. Technical Memorandum No. 899 



19 



although not without difficulties because of the nonuniform 
.processes. Since, as ascertained qualitatively, the ex- 
istence of material speed behind the detonation front is 
possible by the energy law as well as the impulse law in 
case of an extension of the concept regarding the further 
destiny of the combustion gases, equations (21) and (22), 
established for detonation speed and material speed, will 
he considered applicable. 

Assume that figures 5c and 5d are instantaneous photo- 
graphs of the detonation process with a time interval of 1 
second. While within this interval point H of the deto- 
nation front has reached by quantity D the point G, the 
point G of the rarefaction front has traveled the dis- 
tance w to point K. Consequently, the difference of the 
head lengths CH and KG after 1 second, is D - w. The 
increment of the mass of the head is 

— ( D — w ) 

As this mass moves at material speed u, the rise in 
momentum of the head per second is 

-- (E — w ) u 

g 

This rice must equal the effective pressure difference at 
"both ends of the head. This gives the impulse equation 
(p 0 « P a ): 

(34) 

Each second the latent energy p 0 D II is transformed 
into potential and kinetic gas energy of the combustion 
gases. New energy forms per second, are: potential gas 
energy in stage 2, of the order of 

p a w (1 - as p a ) 
k - 1 

the potential gas energy at stage 1, of the order of 

P x (D - w) (1 - a x p x ) 
______ 
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and the kinetic gas energy at stage 1, of the value of 

ig P 1 (D - ») u2 

Hence , 



Po D U 



p x (D ~ w) (1 - a x P x ) 



k - 1 



p 2 « (1 



- a p ) i 



or (equation 14) 



p o 33 f = p i ( D " w) (1 - ct 1 p x ) 
+ P 0 w (1 - a a p 2 ) 



3 

k ~ 1 

2g 



p (D - w) u £ 



(35) 



The law of preservation of mass for an observer sit- 
ting in the rarefaction front gives: 



P 0 w = p (w - u) 

2 X 



(36) 



The solution of the variables w, P 2 , a g , and P g in 
equations (34), (35), and (36) is contingent upon a, Q , 

In the detonation of gaseous explosives, the covol- 
umes ' a x and a 2 are negligible, hence w can be com- 
puted. For gases, it is: 



w = K(k) D 

whereby K(k) is a function of k only: 



(37) 



K(k) 



3 - k + y(3k - l) 2 + 8k (k + l) 
4 (k + 1) 



Putting k roughly at 1.2, we have: 
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K(k) = 0,8 and w = 0.8 D 

j The speed of the rarefaction front in gases is thus 

only about 20 percent lower than the detonation speed or, 
in other words, the head length amounts to about 20 per- 
cent of the distance covered by. the . detonat ion front. 



For solid or liquid explosives whose co volumes 

and a 3 are not ' negl igible , w can "be computed in the 

following manner. As shown "by A. Schmidt (reference 9), 
the covolume of the combustion gases, is a function 

of the gas density p and, consequently, of the original 
density p Q of the explosive. He illustrates this rela- 
tionship as following from detonation-speed measurements 
for different explosive densities with curves and arrives 
at a decrease in covolume for increasing gas density, which 
seems extremely likely. 

Now the covolume is not only dependent on the density 
hut, with constant density, on the temperature as well. 

Visualize a molecule heap of great density hut low 
temperature. Then the individual molecules penetrate into 
their common action spheres merely as a result of their 
need for space (deformation of "molecule envelopes"). If 
a rise in density is accompanied by constant temperature, 
the depth of penetration is greater; that is, the covolume 
becomes less. If, with constant density, the temperature 
is increased, the molecules - as a result of their then- 
increased kinetic energy - penetrate still deeper into the 
envelopes of their neighbors with a further reduction in 
covolume as a result. 

Hence, the factor best suited to describe, as sole 
independent variable, the behavior of the covolume of a 
gas, is the pressure because it is, according to the phase 
equation, proportional to the density as to the temperature, 

The relation a = a (p) obtained, say, for a certain 
gas, may therefore suffice to define the magnitude of a 
when p alone is known. This functional relationship be- 
tween the pressure and covolume of a certain explosive 



a = a 



(p) (38) 



can also be established from detonation-speed measurements 
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for different explosive densities. . The covOlume a x is 
solved according to equation (22), and the correlated 
pressure p x , according to equation (28). Several pairs 

then give the desired curve a = a (p). 

Following this preparatory step, the speed. of the 
rarefaction front is calculable, since the magnitude of 
a s is now defined by the existing pressure p 2 . The four 

available equations, (34), (35), (36 ), and (38) contain 
the four variables w, p g , a 3 , &nd p s . An explicit calcu- 
lation of w would give an expression too complex to han- 
dle. Much quicker results are obtained when computing p g 
according to (36), and p according to (34) for an as- 



sumed w. The a = a(p) curve gives the correlated a value 
for p^. These four variables must satisfy equation (35). 
The process is repeated. 

The temperature T a behind the rarefaction front fol- 
lows from the phase equation at 




(39) 



The quantities in the detonation of trinitrotoluene 
of density 1.59 kg/liter are now computed by means of the 
developed equations, on the basis of Iriedrich's test 
data of detonation speeds at different densities quoted by 
A. Schmidt. 

TABLE I 



Detonation speed of trinitrotoluene in relation to 
density according to Priedrich 



Density 
p Q (kg/liter) 


Detonation speed 
D (m/s) 


1,0 


4,700 


1,29 


5,900 


1.46 


6,500 


1.59 


6,900 
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The heat content is: 

U = 1,085 (cal/kg) 

= 4.63 X 10e (m kg/kg) or (m) 

According to H. Kast (reference 11), the powder con- 
stant f of trinitrotoluene is 

f = 8.08 X 10 4 (m) 

hence the exponent of adiahatic expansion, according to 
equation (14): 

k = 1 + | = 1.175 

(For U the total heat volume, strictly speaking, should 
"be set up to absolute zero point; "but the error caused 
by referring the recorded heat content to 0° C. may he ig- 
nored. ) 

Equation (21) gives the speed of the combustion gases 
(g = 9.81 (m/s a)) 

u = 855 (m/s) 

The factor i or the covolume a, is obtainable 

from equation (23) for the different detonation speeds, 
while the detonation pressure follows from equation (16). 

Table II contains the values thus obtained from 
Priedrichs' measurements. 

TABLE II 

Covolume a x and detonation pressure p x of trinitro- 



toluene according to the developed equations 



Detonation speed 
D (m/s) 




Covolume 
(liters/kg) 


Detonation pressure 
p (kg/ cms ) 


4,700. 




.0.60 . 


4.09X10* 


5,900' 




.53 


6.62X10 4 


6,500 




.49 


8. 27X10* 


6,900 




.46 


9.55X10* 
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The association "between covolume and pressure a = a(p), 
taken from this tahle , is graphically shown in figure 6. 

The explosion pressure is, according to (11 ) 
P E = 7.18 x 10 4 (kg/ cm 2 ) 

Since its solution is contingent upon the knowledge of ovg , 

the calculation must proceed on the "basis of the now-known 
function a = ct(p). The correlated covolume is 

a E = 0.517 (liter / kg) 

The detonation pressure of trinitrotoluene of density 
1.59 kg/liter is thus a"bout 30 percent higher than the ex- 
plosion pressure which, for gaseous explosives, would amount 
to 100 percent. The density of the detonation gases "behind 
the transformation zone follows from equation (31) at 

p x = 1.82 (kg/liter) 

and the explosion temperature from equation (12) at 

or, since the molecular weight M = 32.72, 

T E = 3,120° K. ; T E = 2,847° C. 

The detonation temperature is, according to (33): 

T 1 = 1.08 T E ; T 1 = 3,380° K. ; T 1 = 3,107° C. 

These quantities are compared with Schmidt's findings. 

Except for T the agreement is good. The values 
for the phase quantities "behind the rarefaction front are 
obtained "by approximation. 

Tahle IV gives the quantities involved in the detona- 
tion, along with the explosion quantities for comparison. 
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TABLE III 



Comparison of computed detonation quantities of 
trinitrotoluene of 1.59 kg/liter density and 
6,900 m/ s detonation speed 



Detonation 


Calculated 


according to 




Author *s 






quantity 


Chapman * s 




equations 


equations 




k 


1.175 


1.16 




u 


855 


891 


(m/s) 


a i 


6.46 


0.45 


(liters/kg) 


Pi 


9.55X10 4 


9.65X10 4 


(kg/ cm 3 ) 


Pi 


1.82 


1.83 


(kg/ liter ) 




3,107 


3,630 


(°c.) 



TABLE IT 

Quantities involved in the detonation of trinitrotoluene 
of 1.59 kg/liter density 



Detonation 


Notat ion 


Value 


quan t i t i e s 








Density of solid explosive 


Po 


1.59 


(kg/ liter ) 


Detonation speed- 


B 


6,900 


(m/s) 
(kg/ cm 2 ) 


Explosion pressure 


PE 


7.18X10 4 


Covolume of explosion 






pre s sure 


^E 


0.517 


(liter/kg) 


Explosion temperature 




2 ,847 


(°c.) 


Exponent of adiabatic ex- 






pansion 


k 


1.175 




Material speed of combus- 








tion gases "behind deto- 








nation front 


u 


855 


(m/s) 


Detonation pressure 


P i 


9.55X10 4 


(kg/ cm 3 ) 


Covolume "behind detonation 








front 


a i 


0.46 


(liter/kg) 


Com"bustion-gas density "behind 








detonation front 


Pi 


1.82 


(kg/liter) 

(°c.) 


Detonation temperature 


T i 


3 ,107 


Speed of rarefaction front 


vr 


3 ,720 


(m/s) 


Pressure "behind rarefaction 








front 




5 .03X10 4 


(kg/ cm 2 ) 


(Continued on p. 26.} 
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TABLE IV (Cont. ) 



Quantities involved in the detonation of trinitrotoluene 

of 1.59 kg/liter density 



Detonation 
■quan t i t i e s 


dotation 


Value 


Covolume "behind rarefaction 
front 


a 3 


0.58 (liter/kg) 


G-as density "behind rarefac- 
tion front 


P 2 


1.40 (kg/liter) 


Temperature of gases "behind 
rarefaction front 




2,350 (°C.) 



The ohtained values are illustrated in figure 7. A 
striking feature is the great distance of the two fronts, 
amounting to almost 50 percent of the travel of the deto- 
nation front as compared to about 20 percent for gaseous 
explo give s . 



According to figure 7c, the material speed u pre- 
vails only "between the two fronts, but is zero in the qui- 
escent explosive and during the terminal stage of the com- 
bustion gases. Figure 7d gives the energy density along 
the detonation path, with allowance for the potential as 
well as the kinetic gas energy. The shifting of the energy 
into the detonation head is plainly seen. The integral of 
this energy distribution must naturally equal the total 
latent energy of the explosive released by the detonation 
front F x . Figure 7e illustrates the effect of the mass. 
Here also the integral of the mass distribution must agree 
with the total mass of explosive engendered by the detona- 
t ion. 

Table V gives the individual energies for a column of 
explosive of 1 cm 3 section and D length at detonation. 

The sum of the total energy of the detonation head 
and of the potential energy of the residual combustion 
gases behind the rarefaction front is numerically equal to 
the latent energy. Consequently, the total energy of the 
head is 60 percent of the iatent energy of the explosive 
involved in the detonation. Thereby the kinetic energy 
within the head is considerably smaller than the potential. 
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. TABLE V 






Distribution of energy in a detonating 
trinitrotoluene column of density 1.59 kg/liter 


and length D 


• Energy form 


Magnitude 
(m ke) 


Latent energy 




5.08x10 s 


Potential energy of detonation head 




2.83x10 s 


Kinetic energy of detonation head 




.22X10 8 


Total energy of detonation head 




3 .05X10 8 


Potential energy of gases behind rare- 
faction front 




2 .02X10 8 



DEDUCTIONS FOR THE CONCEPT OF SHATTERING- POWER 



Following proof of the existence of a detonation head 
of increasing length in time as sequence of the physical 
fundamental laws, the question of shattering power of an 
explosive is discussed. 

Visualize at B (fig. 5) a steel plate hit vertical- 
ly by the detonation front P x . 

R. Becker holds the detonation pressure p respon- 
sible for the effect on the steel plate. But a pressure 
alone cannot be decisive for an effect. Assume a closed 
bomb of very small volume under high pressure v 1 • Owing 
to the pressure, the walls expand and the pressure drops 
immediately to a fraction of the initial value, since the 
small energy storage is not able to maintain the initial 
pressure permanently. 

While the pressure alono i's decisive for the effect 
on absolutely rigid walls, the magnitude of all. effects 
exerted on deformable bodies (i.e., all actually existing 
bodies) is contingent upon the available energy. 

Take a closed bomb of large content filled with a gas, 
the pressure of which is lower than the tensile strength 
of the bomb walls, the potential gas energy being suffi- 
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cient to withstand the work of deformation. Even then the 
effect obviously will he zero, on account of the low pres- 
sure . 

The shattering effect therefore does not increase with 
the available energy only hut also with the height of the 
pre s sure . 

The basis of our concept of detonation process now 
suggests consideration of the total energy of the detona- 
tion head as energy reservoir for a shattering effect, 
while the effective pressure is equal to the detonation 
pressure p 1# (The increase in hitting pressure due to 

the combustion-gas momentum is negligible.) 

But the head energy is a function of the amount or 
volume of detonating explosive. Denoting this energy of a 
detonating explosive of 1 cm section and 1 m length with 
Eg, we find for the cited trinitrotoluene 

E K = 4.22 X I0« (i-lK,) 

However, this head energy E K of a certain explosive 
volume cannot as yet be a fitting term for the shattering 
power, "because the detonation pressure itself must, as 
qualitatively deduced, "be present. The higher the detona- 
tion pressure, the higher the detonation speed and the 
smaller the time interval necessary to course through a 1- 
meter column. If this time is t = ^ (s), the shatter- 
ing power B of an explosive appears to be adequately ex- 
pressed by 

B = ~ = Ew- D ) 
t K Sn cm 3 s / 

which presents the shattering power by an energy per vol- 
ume per time; i.e., performance per unit volume. Substi- 
tuting the kilowatt and liter for m kg/ s and cm 2 m, re- 
spectively, we have 

B = E K D (kw/liter) (40) 

Thereby the head energy E K for an explosive of 1-meter 
length and 1-liter volume, i.e., 10 cm 2 section, is com- 
puted in the described manner. 
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The shattering power could equally he referred to the 
unit weight- of the explosive hut since in practice a cer- 
tain volume is available for an explosive, the unit vol- 
ume is more appropriate. 

For trinitrotoluene of 1.59 kg/liter density, it 

give s 

B = 2.98 X 10 7 (kw/liter) 

It is readily seen that the shattering effect of an 
explosive is essentially contingent upon the detonating 
quantity rather than upon the data of detonation pressure 
only, which is not dependent on the quantity of material. 

In the absence of abutment A (fig. 5), the gases in 
stage 2 move under pressure decrease, opposite to that of 
the rarefaction front. The rarefaction advances faster as 
a result of this pressure decrease. 

After a certain detonation path, the distance of both 
fronts therefore is smaller without than with A. Hence, 
since the head energy, despite the constant detonation 
pressure, is dependent on A, it is readily apparent why 
the effects of explosives in the dammed- up stage are great- 
er than those in the free stage. 



CONCLUSION 



In conclusion, we shall analyze the detonation mech- 
anism, The accord of the detonation speed derived from 
the minimum condition with practical experience suggests 
the assumption that the magnitude of detonation speed is 
actually dependent upon the macroscopic quantities only, 
such as heat content, explosion pressure, etc., whereas 
the physical state of the individual explosive particles 
is unimportant, as exemplified by the quality of detona- 
tion speed of solid and liquid explosives of the same den- 
sity. 



The reaction process, which in explosions and combus- 
tions defines the linear speed of transformation, plays no 
part at all in the phenomenon of detonation. 

It may therefore be assumed that the explosive before 
the detonation front is at a stage independent of its ini- 
tial physical state. The author, in consequence, assumes 
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that the radiation of the gages of . comhust ion is sufficient 
to activate a thin layer of the explosive "before the deto- 
nation front and here create the explosion pressure p B 
and the explosion temperature T E . ^Then the magnitude of 
the detonation speed is given only "by the propagation 
speed of a shock wave in a gas of pressure p E and tempor- 



A molecular-kinetic analysis of the detonation phenom- 
enon (reference 12) can he made "backward from the already 
"available thermodynamic data precisely as the nature of the 
velocity of sound on the "oasis of the known thermodynamic 
formula was recognized as the quantity of the molecular 
speed of the particular gas multiplied hy a certain con- 
stant . 

Translation "by J. 7anier, 
National Advisory Committee 
for Aeronautics. 
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Figure 1.- Detonation process. 
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Figure 3.- Improved detonation 
process. 
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Figure <t.- Detonation speed plotted 
against speed of 
combustion gases. 
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Figure 5.- Formation of 
detonation 

head. 
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Figure 6.- Covolume of 

trinitrotoluene 
plotted against pressure. 
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Figure V.- Distribution 
of pressure, 
material speed, energy 
density, and tempera- 
ture in a detonating 
column of trinitroto- 
luene of D length. 
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